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Wave  Structures  in  Thermospheric  Density 
from  Satellite  Electrostatic  Triaxial 
Accelerometer  Measurements 


1.  INTRODUCTION 


The  frequent,  occurrence  of  wavelike  perturbations  in  thermospheric  and 

ionospheric  constituent  densities,  winds,  and  temperatures  is  now  universally 

accepted,  although  the  origin  and  propagation  characteristics  of  these  oscillations 

12 

are  poorly  understood  and  still  under  investigation.  *  Historically,  the  first 

3  4 

comprehensive  studies  ’  involved  the  analysis  of  ionospheric  data,  and  the  wave 
motions  were  referred  to  as  'travelling  ionospheric  disturbances'  (TIDs)  although 
it  was  known  that  fluctuations  in  ionospheric  parameters  were  actually  driven  by 

3 

collisions  with  the  oscillating  neutral  gas,  Georges  categorized  these  waves 


(Received  for  publication  for  publication  2  June  1987) 
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according  to  their  horizontal  wavelengths  as  being  either  large-scale  (1000  to 
4000  km)  or  medium-scale  ( 100  to  500  km).  The  large-scale  waves  are  also  charac¬ 
terized  by  periods  of  order  1  to  3  hours  and  propagation  velocities  of  order  0.  5  to 
1.0  km/sec,  and  are  thought  to  be  primarily  excited  during  substorms  by  im¬ 
pulsively  enhanced  particle  precipitation  heating,  Lorentz  forces,  and  joule  dissi¬ 
pation  at  high  latitudes.  ’  ’  The  large-scale  waves  can  propagate  long  distances 

8  9 

from  the  source  and  are  detectable  at  low  latitudes.  ’  The  medium-scale  waves 

propagate  slowly  (100  to  250  m/sec)  and  have  periods  in  the  range  of  15  to  60  minutes. 

Their  origin  is  less  certain,  but  they  have  been  tenuously  connected  with  meteoro- 
■  .  ,  ,  10,11,  12,  13 

logical  sources  in  the  troposphere.  Medium  -scale  waves  propagating 

horizontally  away  from  an  auroral  source  do  not  travel  long  distances;  however,  it 

is  thought  that  medium-scale  waves  partially  reflected  at  the  base  of  the  thermos- 
14  25 

phere  or  totally  reflected  at  the  surface  may  be  ducted  into  the  thermosphere 

1  2  16  17 

with  little  attenuation  and  travel  large  distances  from  the  source.  ’  ’  ’ 
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The  experimental  view  based  on  wave  analyses  of  neutral  and  ion  density  and 
temperature  measurements  from  the  AE-C  satellite  is  that  a  quasi -stationary, 
quasi -persistent  wave  structure  (scales  of  order  400  to  4000  km)  that  mav  slowly 
dissipate  and  be  regenerated  by  auroral-region  sources  exists  over  global  distances 
in  the  thermosphere. 

The  above  theoretical  studies  and  e  cpe rimental  results  suggest  a  need  for  wave 

analyses  of  satellite  measurements  that  focus  on  the  latitudinal  and  magnetic  ac- 

19 

tivity  dependences  of  the  wave  characteristics.  An  existing  data  base  consisting 
o1’  low-altitude  (160  to  220  km)  densities  measured  from  a  series  of  Air  F’orce 
satellite  accelerometer  experiments  is  ideally  suited  for  this  purpose,  as  there 
exist  many  consecutive  orbits  of  data  ranging  from  high  to  low  latitudes  in  a  high- 
inclination  orbit  ('  80  to  85M  covering  both  quiet  and  active  magnetic  conditions 
between  1979  and  1985.  Density  data  from  the  Satellite  Electrostatic  Triaxial 
Accelerometer  (SETA)  Experiment  were  spectrally  analyzed  by  the  Maximum 
Entropy  Method  (MEM).  The  purpose  is  to  make  available  a  technique  by  which  the 
occurrence  of  wavelike  structures  in  thermospheric'  densities  can  be  quantitatively 
analyzed  and  correlated  with  various  geophysical  conditions,  including  the  deposi¬ 
tion  of  energy  at  high  latitudes  during  geomagnetically  disturbed  conditions.  Spectra 
were  computed  for  tvpical  magnetically  quiet  and  active  days  during  July  1983,  and 
the  relative  occurrence  of  spectral  peaks  in  wavelength  ranges  characteristic  of 
large-scale  and  medium-scale  TIDs  is  examined  as  a  function  of  local  time  and 
latitude. 

2.  DESCRIPTION  OF  DATA  AND  PROCESSING  STEPS 
2.1  The  SETA  Experiment 

The  Satellite  Electrostatic  Triaxial  Accelerometer  Experiment  (SETA)  is  based 
on  the  single-axis  MESA  instrument  flown  by  the  Air  Force  Geophysics  Laboratory 
on  eight  satellites  since  1968  including  Atmosphere  Explorers  C,  D,  and  E.  The 
MESA  utilizes  a  proof  mass  in  the  shape  of  a  flanged,  hollow  cylinder.  This  proof 
mass  is  suspended  electrostatically  along  its  radial  axes  and  is  force -rebalanced 
electrostatically  along  its  longitudinil  or  sensitive  axis.  A  complete  description 

18.  Gross,  S.  H.  ,  Reber,  C.  A.  ,  and  Huang,  F.  T.  (1984)  La rge -scale  waves  in 

the  thermosphere  observed  bv  the  AE-C  satellite,  IEEE  Trans,  on  Geos, 
and  Remote  Sens.,  GE -22;340. 

19.  Marcos,  F.  A.  (  1987)  Requirements  for  improved  thermospheric  models, 

J.  Astronaut.  Sci.,  in  press. 


of  the  MESA  ms‘rument  flown  on  the  three  Atmosphere  Explorer  satellites  is  uiveti 

20 


hv  Champion  and  Marcos."  To  develop  a  triaxial -sensing  caoahilitv,  'hi'  susnen 
sion  system  was  modified  by  instrumenting  both  radial  axes  w;'h  p-eeision  eon  - 


strainmeo*  loops.  A  complete  description  of  the  SETA  configuration  is  given  hv 


2  1  21 
Marcos  and  Swift.  “  Marcos  and  Swift”  also  describe  results  >f  an  orbital  High’ 


test  from  which  various  instrument  response  characteristics  and  design  improve¬ 
ments  were  obtained. 

The  accelerometer  output  contains  contributions  related  to  vehicle  dynamics 
and  instrument  operating  characteristics  as  well  as  those  due  to  aerodynamic  drag 
and  winds.  Cen  ripetal  acceleration  and  the  gravity  gradient  effect,  each  of  the 
order  of  0.22  pg  but  oppositely  directed,  provide  a  negligible  0.01  pgt  resultant 
dynamic  acceleration  in  the  along-track  (Z)  ax's.  Extraneous  accelerometer  out¬ 


puts  consist  of  a  bias  a^  and  random  noise  a^..  With  the  centripetal  and  gravity 


gradient  accelerations  incorporated  as  part  of  the  bias  term  and  the  vehicle  dynamic 


acceleration  represented  by  a^  the  average  accelerometer  output  over  the 


2.045-sec  sample  time  a^.  is  then  represented  by 


a„. .  =  a . .  t-  a  . 
It  1)1  V’ 


f-  a 


13  i 


Ni 


(1) 


where  i  =  X,  Y,  Z,  since  the  various  acceleration  sources  can  be  different  along 
each  axis.  The  aerodynamic  component  is  derived  using  Eq.  (1)  in  the  form 


al)i  =  aTi 


(a  . 
vi 


aBi  +  aNi)  * 
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Numerical  filtering  techniques  are  used  to  remove  frequency  components  of 
a.j.  due  to  satellite  dynamics  a^  and  instrument  noise  a^..  The  cutoff  frequency  of 
th*  tow-pass  filter  is  0.066  Hz  (see  Section  2.2).  Since  the  along-track  accelerometer 
axis  senses  aerodynamic  drag  throughout  the  orbit,  bias  cannot  be  directly 
deduced  by  the  accelerometer  output  at  maximum  altitude.  A  bias  value  is  estimated 


20.  Champion,  K.S.  W.  ,  and  Marcos,  K.  A.  (1073)  The  triaxial  accelerometer 
system  on  Atmosphere  Explorer,  Radio  Sci.,  8:197. 

2  1.  Marcos,  K.  A.  ,  and  Swift,  E.R.  (1082)  Application  of  the  Satellite  Triaxial 
Accelerometer  Experiment  to  Atmospheric  Density  and  Wind  Studies, 
AEGL-TR-62  -000  1,  AI)  A  120852.  1 

22.  Noonan,  J.P.  ,  Eioretti,  K.W.  ,  and  Hass,  13.  (  1075)  Digital  Filtering  Analysis 
Applied  to  the  Atmosphere  Explorer-C  Satellite  MESA  Accelerometer  Da'a, 


for  each  orbit  essentially  by  subtracting  a  drag  acceleration  based  on  a  model 
23 

density  value  from  the  filtered  accelerometer  output  at  maximum  altitude.  This 
density  typically  corresponds  to  a  drag  value  of  0.4  pg,  whereas  the  drag  at 
minimum  altitude  (~  170  km)  is  typically  of  the  order  of  8  pg.  For  the  cross-track 
and  radial  axes,  aerodynamic  accelerations  are  negligible  at  maximum  altitude. 

The  equations  used  for  the  determination  of  atmospheric  density  and  wind 
velocities  are 


C  p  V2 
2  m  t K 


where 


axis  X,  Y,  Z; 

aerodynamic  acceleration  along  axis  i; 

satellite  frontal  area; 

satellite  mass; 
atmospheric  density; 

atmospheric  mass  velocity  relative  to  the  satellite. 


V  =  -VG  +  V  A  +  VW 


where 


inertial  satellite  velocity; 

atmospheric  rotation  velocity  (assumed  equal  to  earth's 
rotation  velocity); 
neutral  wind  velocity; 

drag  coefficient  for  axis  i;  C.  =  CdV^,  Vy ,  VL,), 

24 

calculated  from  the  formulation  of  Sentman.  “ 


23.  Jacchia,  L.G.  (1971)  Revised  Static  Models  of  the  Thermosphere  and 


Exosphere  With  Empirical  Temperature  Profiles,  Spec.  Rept. 


Smithsonian  Astrophys.  Observatory,  Cambridge,  Massachusetts. 

Sentman,  L.  H.  (1961)  Free  molecular  flow  theory  and  its  application  to  the 
determination  of  aerodynamic  forces,  Rept.  EMSC-4485  14,  Lockheed 
Missiles  and  Space  Co.  ,  Sunnvvale,  California. 


In  Eq.  (2)  only  three  acceleration  measurements  are  available  (a.),  while  there 
are  four  unknowns  (p,  V).  An  iterative  technique  has  been  developed  to  determine 
densitv  and  winds.  The  known  contributions  to  the  Z  axis  velocity  \  (that  is, 

f  V  are  large  compared  to  the  uncertaintv  \’w^;  hence  an  initial  estimate 
of  the  densitv  is  obtained  by  assuming  Vq=  -Y^,  +  V  ^  and  calculating 


1  Aref 


Vcz<vo> 


This  density  estimate  is  used  to  calculate  the  three  velocity  components  (V) 

from  Eq.  (2)  and  to  obtain  the  wind  as  V'w  =  V  -  VQ.  To  refine  the  solution,  a 

new  V  =  (V,..  V  ,  V  )  is  used  in  Eq.  (4)  to  calculate  a  new  p.  This  new  p  is  then 
.A  Y  U  ca 

used  to  calculate  new  velocities,  and  so  on,  until  convergence  is  achieved.  The 
wind  velocity  and  density  estimate  are  accepted  as  final  if  |v  <  0.002  \\' y  I. 

The  errors  in  the  calculated  values  of  density  are  attributed  to  accele rometer 
bias,  attitude  uncertainties,  drag  coefficient  estimates,  noise  within  the  desired 
signal  frequency  spectrum,  and  the  inseparability  of  density  and  winds  along  the 
in -track  axis.  Systematic  errors  due  to  the  drag  coefficient  are  estimated  to  be 
within  ±  5  percent,  while  those  due  to  bias  varv  from  ±  0.  4  percent  at  170  km  to 
±2  percent  at  200  km.  Low -frequency  (<0.005  Hz)  noise  of  amplitude  0.04  pg 
contributes  a  random  error  of  ±0.5  percent  at  170  km  and  ±  2.  5  percent  at  200  km. 

In  addition,  an  assumed  wind  uncertaintv  of  200  m/sec  along  the  Z  axis  contributes 
a  random  error  in  p  of  .-5  percent.  Attitude  errors  for  all  three  axes  are  negligible. 


222  Pre-Processing  (Sampling,  Filtering,  and  so  on) 

The  SETA  measurements  were  digitized  at  the  rate  of  approximately  one  point 
per  2  seconds.  The  data  used  in  the  present  analysis  were  taken  at  a  rate  of  one 
point  per  4  sec,  corresponding  to  a  Nvquist  frequency  of  0.  125  Hz.  A  low- 
pass  filter  was  applied  to  the  data  with  a  cutoff  frequency  of  1/  15  Hz  and  unit 
response  at  1/20  Hz.  The  maximum  identifiable  frequency  in  the  data  series  there¬ 
fore  occurs  in  the  range  0.05  to  0.07  Hz,  Since  the  satellite's  velocity  is 
7.8  km/sec,  this  translates  to  a  minimum  horizontal  wavelength  of  order  100  to 
150  km.  Note  that  since  the  group  velocities  of  Tin's  are  generally  less  than 
0.7  5  km  (sec,  the  measurements  reported  here  can  be  viewed  as  representative 
of  the  near -instantaneous  horizontal  wave  structure  within  a  given  latitude  segment. 
H  >w  -vt-r,  since  the  satellite  orbital  period  is  00  min,  no  information  can  be 
deduced  regarding  the  waves'  frequencies  or  periods.  When  the  term  'wave 


frequency'  is  used  in  this  report,  one  is  actually  referring  to  an  effective  frequency 
implied  by  the  movement  of  the  satellite  through  a  fixed  spatial  wave  structure. 

The  SETA  experiment  covers  the  170  to  220  km  height  region  between  about 
0°  and  84"  geographic  latitude  in  a  morning/ evening  (~  1030  LT / 2  130  LT)  orbit. 
Before  applying  the  MEM  (see  Section  2.3),  the  downleg  (evening)  and  upleg 
(morning)  portions  of  the  orbit  are  divided  into  the  following  segments: 


Approximate 

Local  Time 

Latitudes 

2130 

10° 

-  35° 

2130 

30° 

-  55° 

2130 

55° 

O 

O 

CO 

I 

1030 

0 

O 

CO 

-  55" 

1030 

55° 

-  30° 

1030 

35° 

-  10" 

This  was  done  to  allow  examination  of  differences  in  thermospheric  wave  structure 
as  a  function  of  latitude  and  Local  time,  particularly  during  geomagnetic  disturb¬ 
ances.  The  25°  latitude  segments  correspond  to  slightly  less  than  400  sec  of 
data,  or  approximately  3000  km  (or  about  0.0025  Hz  in  frequency  space).  A  high- 
pass  filter  was  not  applied  to  the  data.  The  measurements  were  normalized  and 
de-trended  (using  a  linear  least-squares  fit)  within  each  latitude  segment  before 
applying  the  MEM. 

23  Spectral  Analysis  by  the  Maximum  Entropy  Method  (MEM) 
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Since  1975,  AFGL  has  been  applying  MEM  to  Air  Force  problems.  In  each 
case,  for  each  different  type  of  data  set,  MEM  has  been  compared  to  the  "classical1 
methods,  either  Blackman- Tukey  [BT]  or  the  Cooley-Tukey  periodogram  method 
using  the  Discrete  Fourier  Transform  [DFT|.  In  every  case,  MEM  has  been  shown 
to  be  demonstrably  superior  to  the  classical  methods.  The  superiority  lies  in  both 
increased  resolution  and  increased  smoothness  in  the  frequency  domain  and,  for 
slowly  varying,  nearly  stationary  signals,  increased  resolution  in  the  time  domain. 
For  slowly  varying,  non -stationary  signals,  a  time  window  can  be  moved  through 
the  data  set  in  such  a  way  that  each  snapshot  portrays  a  very  nearly 

25.  Radoski,  H.  R.  ,  Fougere,  P.  F.  ,  and  Zawaliek,  E.J.  (  1975)  A  comparison  of 
power  spectral  estimates  and  applications  of  the  maximum  entropy  method 
J.  Geophys.  Res.,  80:619. 


The  order  m  of  the  MEM  determines  the  smoothness  or  c  »mplexity  of  the 

spectrum.  For  small  m  the  spectrum  vviLl  be  smooth  and  simple;  as  m  increases, 

the  resolution  and  complexity  of  the  spectrum  will  increase.  The  Akaike.  or 

27 

final  prediction  error  (FPE)  criterion  for  selecting  the  MEM  order,  requires 

that  m  which  minimizes  the  FPE.  However,  it  has  been  found  in  the  satellite 

1 3 

density  data  analyzed  bv  Gross  that  the  FPE  exhibits  a  broad  and  shallow 
minimum  in  all  cases.  For  our  purposes,  we  wish  to  have  the  simplest  spectrum 
that  identifies  all  the  major  spectral  peaks  characteristic  of  the  data.  The  order 
was  therefore  chosen  by  visually  examining  data  where  certain  frequencies  were 
obviouslv  exhibited,  and  bv  increasing  m  until  the  expected  spectral  peaks  were 
adequately  identified;  further  increasing  m  beyond  this  point  did  not  change  the 
dominant  spectral  peaks,  but  only  added  structure  to  the  spectrum.  An 
example  of  varying  the  order  m  on  the  spectra  for  a  given  segment  of 
data  is  illustrated  in  Figure  3  for  Rev  54  1  on  Day  205,  a  period  of  high  magnetic 
activity.  We  varied  the  order  m  from  5  to  30  in  steps  of  5.  On  the  basis  of  these 
results,  and  for  several  other  typical  orbits  which  we  have  examined,  a  value  of 
m  =  25  was  chosen  for  all  the  spectra  to  be  presented  here. 


27.  Ulryck,  T.J.  ,  and  Bishop,  T.  N.  (1975)  Maximum  entropy  spectral 

analysis  and  autoregressive  decomposition.  Rev.  Geophys.  Space  Phys.  , 
13:183.  - - 


3.  ANALYSIS  AND  INTERPRETATION  OF  SPECTRA 


In  this  section,  spectra  computed  for  typical  magnetically  quiet  and  active  days 
are  presented,  and  their  characteristics  analyzed  as  a  function  of  local  time  and 
latitude.  Specifically,  the  data  are  taken  during  Days  203  and  205  of  1983,  which 
occurred  during  July.  The  3 -hourly  Kp  indices  are  shown  in  the  table  below; 


Day  Kp 


203 

2  + 

lo 

lo 

lo 

1- 

1- 

3- 

5- 

205 

6- 

5o 

4o 

4- 

6o 

4o 

4o 

3  + 

Sample  normalized  density  data  plots  for  Revs  511  (during  the  first  3  hours  of 
Day  203)  and  542  (during  the  first  three  hours  of  Day  205)  are  shown  in  Figures 
4  through  7.  The  data  are  separated  according  to  day  (Figures  5  and  7)  and  night 
(Figures  4  and  6),  and  according  to  latitude  regions  (10°  to  35°,  30°  to  55’,  and 
55°  to  80°).  The  ranges  of  altitudes  are  also  shown  on  the  plots.  The  Rev  542 
data  generally  exhibit  an  enhanced  degree  of  structure  and  waviness,  presumably 
connected  with  the  increase  in  magnetic  activity.  The  corresponding  maximum 
entropy  spectra  are  shown  in  Figures  8  through  11.  The  tick  marks  on  the  spectral 
plots  represent  the  40-dB  points  for  the  power  spectral  density.  Note  that  while 
the  spectral  power  is  considerably  enhanced  during  the  magnetically  active  period, 
and  some  frequency  shifting  of  the  peaks  occurs,  the  structure  of  the  spectrum 
does  not  appear  to  be  appreciably  enhanced  as  a  result  of  the  increased  geomagnetic 
activity. 

The  relationship  between  wave  spectral  characteristics  and  geomagnetic 
activity,  latitude,  and  local  time  requires  examination  of  a  much  larger  data  base, 
however,  including  examination  of  seasonal  effects  as  well.  A  comprehensive 
statistical  and  scientific  analysis  is  beyond  the  scope  of  this  report.  However,  we 
do  analyze  the  first  14  revs  (2  1  hours)  of  Days  203  and  205  of  1983,  and  demonstrate 
the  type  of  study  which  can  be  performed  with  the  data  analysis  procedure  which 
has  been  set  up.  Specifically,  we  will  examine  Revs  509-522  of  Day  203  and  Revs 
541-554  of  Day  205  for  the  relative  occurrence  of  various  scale  waves  whose- 
spectral  density  peaks  exceed  the  30-dB  and  40-dB  levels  as  a  function  of  local 
time,  latitude,  and  magnetic  activity.  Tabulations  representing  percent  occurrence 
of  peaks  exceeding  30  dB  and  40  dB  are  listed  respectively,  in  Tables  I  and  2. 


(Note  that  since  occasional  data  gaps  occurred,  we  have  listed  percent  occurrence, 
rather  than  number  of  occurrences;  however,  since  data  gaps  were  fairly  rare, 

100  percent  occurrence  corresponds  to  no  less  than  12  occurrences  during  the  day.  ) 

Examination  of  Tables  1  and  2  reveals  the  following.  First,  from  Table  2, 
we  see  that  the  most  commonly  occurring  waves  are  those  characterized  by  large 
horizontal  scales  (1250  to  2500  km).  There  also  exists  evidence  of  an  increase  in 
occurrence  for  enhanced  magnetic  activity  conditions  [with  the  exception  of  low 
latitudes  (10°  to  35°)  during  the  morning  (1030  h)].  During  the  quiet  day  the  only 
occurrences  of  medium-scale  waves  were  at  high  latitudes  (55s  to  80°)  during 
nighttime  (2030  h).  This  pattern  persists  during  the  magnetically  active  day  with 
about  twice  the  probability  (15  to  23  percent  as  opposed  to  8  to  15  percent).  En¬ 
hanced  wave  activity  at  high  latitudes  during  nighttime  is,  in  fact,  consistent  with 

28  29 

model  simulations  and  data  for  the  CDAW-6  magnetic  storm  period.  ’  Further, 
during  the  magnetically  disturbed  day  the  data  in  Table  2  indicate  an  increased 
preference  for  the  existence  of  large-scale  waves  as  one  progresses  from  low  to 
high  latitudes. 

The  percent  occurrences  oc  spectral  density  peaks  in  excess  of  30  dB  are 
listed  in  Table  1.  Also  listed  at  the  bottom  of  the  table  are  differences  in  percent 
occurrence  (Day  205  minus  Day  203).  In  almost  all  cases,  waves  of  both  "large" 
and  "medium"  scale  occur  more  often  than  20  percent  of  the  time.  As  a  result 
of  increased  magnetic  activity,  a  general  enhancement  o1'  wave  activity  occurs  for 
large  (  1250  to  2500  km)  and  medium-scale  (150  to  500  km)  waves.  We  speculate 
that  the  increased  occurrence  of  medium-scale  wave  activity  at  low  and  middle 
latitudes  during  magnetically  disturbed  periods  may  be  attributable  to  ducted  waves 
in  the  low'er  thermosphere  as  suggested  by  the  theoretical  calculations  of  Mayr 
et  al.  lf  2 


28.  Roble,  R.G.,  Forbes,  .I.M.,  and  Marcos,  F.  A.  ( 1 987 )  Thermospheric 

dynamics  during  the  March  22,  1  >79  magnetic  storm  (1)  Model  simulations, 
J.  Geophys.  Res.,  J92:60'15. 

29.  Forbes,  .T.M.,  Roble,  R.G.,  and  Marcos,  F.  A.  ( 1987 )  Thermospheric 

dynamics  during  the  March  22,  197!*  magnetic  storm  (2.  )  Comparisons  of  m 
modei  predictions  with  observations,  .1.  Gcophvs.  Res.,  92:6069. 
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Figure  11.  MEM  Spectra  (power  spectral  density  vs  frequency)  of 
Order  =  25  for  Daytime  (~  1030  h)  Portion  of  Rev  542.  Tick 
marks  correspond  to  40 -dB  points  for  power  spectral  density 


Table  1.  Percent  Occurrence  of  =  30  dR  Spectral  Peaks  as  a  Function 
of  Horizontal  Wavelength 


1 - — - - - 1 

Latitude 

Local 

T  ime 

"Large 

Horizontal  Wavele 

-Scale" 

ngth  (km ) 

"’Vied  um  S 

'ale” 

1250- 

2500 

750- 

1250 

500- 

750 

250- 

500 

150- 

250 

(DAY  203) 

10°  -  35’ 

2030 

46 

31 

23 

69 

54 

30 5  -  55’ 

2030 

62 

23 

3  1 

46 

23 

55’  -  80’ 

2030 

7  9 

14 

43 

7  1 

29 

10’  -  35’ 

1030 

92 

15 

39 

15 

15 

30’  -  55’ 

1030 

64 

36 

36 

3  6 

•) 

55’  -  80’ 

1030 

77 

15 

62 

23 

0 

(DAY  205) 

10’  -  35’ 

2030 

67 

17 

42 

67 

92 

30’  -  55’ 

2030 

77 

23 

46 

85 

35 

55’  -  80’ 

2030 

92 

8 

62 

77 

77 

10’  -  35’ 

1030 

62 

23 

38 

54 

38 

30’  -55’ 

1030 

83 

33 

42 

92 

42 

55’  -  80’ 

1030 

100 

18 

27 

73 

36 

(Differences, 

DAY  205-203) 

10°  -  35’ 

2030 

+  2  1 

-14 

+  19 

-2 

+  38 

30’  -  55’ 

2030 

H5 

0 

+  15 

+  39 

+  62 

55’  -  80’ 

2030 

M3 

-6 

+  19 

+  6 

+  48 

10’  -35’ 

1030 

-30 

+  8 

-  1 

+  49 

+  23 

30’  -55° 

1030 

+  19 

-3 

+  6 

+  56 

+  33 

55’  -80’ 

1030 

+  23 

+  3 

-35 

+  50 

+  36 

Table  2.  Percent  Occurrence  of  5  40  dB  Spc  -tral  Peaks  as  a  Function 
of  Horizontal  Wavelength 


"Large 

Horizontal  Wavelength  (km) 

Scale"  "Medium 

Scale" 

Latitude 

Local 

Time 

1250- 

2500 

750- 

1250 

500- 
7  50 

250- 

500 

150- 

125 

(I)AY  203) 

10’  -  3 r> ° 

2030 

1 5 

0 

0 

0 

0 

30’  -  55’ 

2030 

3  1 

8 

0 

0 

0 

55’  -  80 ’ 

2030 

46 

0 

15 

8 

3 

10’  -  35’ 

1030 

3  8 

0 

0 

5 

0 

30°  -  55 J 

1030 

8 

0 

0 

0 

0 

55’  -  80° 

1030 

15 

1  > 

0 

0 

0 

(DAY  205) 

10’  -  35° 

2030 

25 

0 

8 

0 

•» 

) 

30'  -  55  ’ 

2030 

46 

15 

0 

!! 

0 

55’  -  80' 

2030 

6!) 

8 

23 

15 

15 

10’  -  35’ 

1030 

23 

3 

0 

0 

0 

30’  -  55’ 

1030 

38 

3 

0 

0 

0 

55’  -  80’ 

1030 

54 

0 

8 

0 

0 

4.  CONCLUSIONS 


Doisity  data  from  the  Satellite  Electrostatic  Triaxial  Accelerometer  (SETA) 
experiment  were  spectrally  analyzed  by  the  Maximum  En'ropy  Method  (MEM). 
Spectra  jvero  computed  for  typical  magnetically  (quiet  and  active  days  during  duly 
1083,  and  the  relative  occurrence  of  spectral  peaks  in  wavelength  ranges  of  large- 
scale  (hundreds  to  thousands  of  kilometers)  TIDs  is  examined  as  a  function  of  local 
time  and  latitude.  We  preliminarily  conclude  from  analysis  of  these  data  the 
following: 

(a)  Wave  motions  contribute  to  the  observed  variability  of  density 
during  both  quiet  and  active  conditions,  and  become  increasingly 
important  as  magnetic  disturbances  become  more  intense, 

(b)  There  exists  evidence  for  enhanced  wave  activity  at  high 
latitudes  and  during  nighttime. 


V*  *T7Y»VW>V* 
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(c)  There  exists  a  significantly  higher  probability  of  finding  waves 

of  horizontal  wavelengths  between  1250  to  2000  km  and 

150  to  500  km,  especially  during  geomagnetically  disturbed 

conditions,  as  opposed  to  those  waves  with  wavelengths  between 

500  and  1250  km.  The  former,  larger-scale  waves,  are  known 

to  be  capable  of  propagating  large  distances  with  little  attenuation. 

We  tentatively  conclude  that  the  relatively  frequent  occurrence 

of  the  shorter-scale  waves,  even  at  low  latitudes,  may  be  attributable 

to  ducting  in  the  lower  thermosphere  which  allows  them  to  travel 

1  2 

large  distances  as  suggested  by  Mayr  et  al.  ’ 

These  preliminary  findings  can  be  verified  and  further  refined  through  analysis 
of  additional  satellite  density  data  now  in  existence.  The  present  report  provides 
a  computerized  procedure  to  process  these  data  to  reveal  the  dominant  wave 
characteristics,  and  a  rationale  and  direction  to  future  scientific  efforts  along 
these  lines. 
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